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Abstract 

We introduce a simple model of active transport for an ensemble of particles driven by an external 
shear flow. Active refers to the fact that the flow of the particles is modified by the distribution 
of particles itself. The model consists in that the effective velocity of every particle is given by the 
average of the external flow velocities felt by the particles located at a distance less than a typical 
radius, R. Numerical analysis reveals the existence of a transition to clustering depending on the 
parameters of the external flow and on R. A continuum description in terms of the number density 
of particles is derived, and a linear stability analysis of the density equation is performed in order 
to characterize the transitions observed in the model of interacting particles. 

PACS numbers: 05.45.-a, 05.60.-k 
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I. INTRODUCTION 



Two different types of transport problems can be roughly distinguished: passive and ac- 
tive. The case of passive transport occurs when the transported quantity does not affect 
the advecting flow, as exemplified by a dye immersed in a fluid, or of any reacting sub- 
stance like a chemical pollutant having no feedback on the carrying atmospheric or oceanic 
flow 1]. Conversely, in the active transport problem, the subject of this paper, the flow 
itself is modified by the advected substance. Sometimes this is also refered as self-consistent 
transport since the velocity field is in general determined by the substance via a dynamical 
constraint [^. The temperature field, an ensemble of charged particles moving in a self- 
generated electric field, the vorticity of a fluid flow, and gravitationally interacting particles, 
are a few examples reflecting the ubiquity and relevance of active transport processes in 
Nature. 

Recently, much progress has been achieved in both self-consistent and passive processes 
through its reformulation within the Lagrangian description 0,0], which studies transport 
in terms of individual particle trajectories instead of scalar fields. Thus, the Lagrangian 
description of a non-reacting passive scalar in an external velocity field, v(x, t), is given by 

^ = v(x,t) + v^r/(t), (1) 

where Do is the diffusion coefficient of the passive scalar, and T7 is a normalized Gaussian 
white noise with zero mean and delta correlated in time. In eq. the passive character 
is shown in the fact that there is no coupling between the equations of motion. On the 
contrary, in active transport the interactions among particles alters the trajectory of any 
of them^o that for an ensemble of particles immersed in a fluid flow one can write in 
general p| 

^ = v(xi(t),...,x^(t)), (2) 

i = 1, ...,N. This A^-body problem is often treated in a mean-field approximation where 
every particle is considered independently of the rest but in an average potential determined 
self-consistently from the motion of all the particles ^ |^ . In this way, the influence of any 
particle on the system is just through its contribution to the potential. 

In this work, we introduce a distinct type of self-consistent transport model. At difference 
of the mean-field approach, our model assumes a finite range of interaction, R, for any 
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particle, so that particles only interact with others surrounding them. Eq. Q takes the 
form 

= v(x,(t),Xi+i(t)...,Xi+iVB(i)(t)), (3) 

where iVij(z) is the number of particles at distance less than R of particle i, and with 
i + + Nji{i) we label these particles. Most importantly, the external flow is given 

and the invidual particles modify their response to the flow according to the local density 
around them. Like the mean-field one, our model is an intermediate case between the passive 
transport case and the many-body self-consistent models with real interactions decaying with 
distance. The aim of the present paper is to show that even very simple active systems can 
show a very rich behaviour, and, in particular, the formation of clusters of particles may 
appear. Also, due its simplicity one can present a detailed analytical study of the model, 
and show that clustering emerges as a deterministic instability of the density equation of 
the system. 

The paper is organized as follows. In the next section we introduce the model and present 
numerical results showing the clustering. Then, in Sec. Ill we derive the density equation 
for the dynamics of the particles, and perform a linear stability analysis of this continuum 
description. Then, we finish in Sec. IV with the summary of the work. 



II. SELF-CONSISTENT MODEL OF PARTICLES DRIVEN BY AN EXTERNAL 
SHEAR FLOW. NUMERICAL RESULTS 

Let us consider particles in a two-dimensional system of size L x L, and the pres- 
ence of a stationary incompressible two-dimensional shear flow \{x,y) = (0,f(x)). In the 
model the effective velocity of particle i at time t, w^^^ {t), which is in the position Xj(t), is 
the average velocity of the external velocities felt by the particles in its /^-neighbourhood. 
Mathematically 



d'X.i{t) _ ^ e// 
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<"it), (5) 



dt 

where, as indicated above, Nji{i) denotes the number of particles at a distance less than R of 
particle i, and the sum is restricted to the particles j such that \^i(t) —yij(t)\ < R. Periodic 
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boundary conditions are considered, and finite-size effects of the particles, like inertia and 
collisions, are neglected. Note that the self-consistent character of the model comes from the 
fact that at every time the velocity of any particle is determined by the (local) distribution 
of particles itself. A noise term similar to eq. ([Q) could be added to the r.h.s of eq. (j^, 
but this is not considered in this work and we just suppose that advection induced by the 
external flow dominates on the random motion of the individual particles. 

Two limits are clearly identified, R —>■ is the tracer limit, i.e., every particle is simply 
driven by the flow. In the opposite i? — > L all the particles move with the same velocity, 
which is just an average of the external velocity field over all the particles in the system. 
Physically, the model mimicks particles transported by a flow and with some kind of effec- 
tive non-local interaction that force them to move locally with the same velocity. In the 
context of living organisms, traffic or behaviour of human societies many different models 
have been proposed where the density of particles modify their velocity Q]: repulsion, at- 
traction, distribution of resources, cooperation, are some of the types of interactions among 
the individuals that are usually studied. These interactions are mediated (in a biological 
framework) through vision, hearing, smelling or other kinds of sensing, which is reflected, 
as in our model, by the appearence of a typical interaction radius, R. However, a crucial 
difference of these biologically oriented models with (^HS)) is that in those the particles are 
self-propelled, i.e., they have their own velocity. In our model, the velocity is externally 
given, and it is our aim to study the properties of the system of particles depending on 
the characteristics of the external flow. Regarding a biological motivation, our model is 
adequate for acuatic organisms moving by the water flow that modify their velocity as a 
response to other individuals living within a certain distance. 

Concerning the clustering properties of the model, which is the main focus of this work, 
it is clear that shear enhances encounters among particles, and this favours that particles 
group together due to the averaging of velocities. On the contrary, local strength of the 
external velocity field tends to disperse the particles, breaking clusters. Combining these 
two effects a typical length scale is introduced: 

^ < {dv{x)/dxf > . 

< v{xy > ' 

Here < . >= 1/L dx., and A is related to the Taylor microscale of turbulence, though 
here the meaning is somewhat different since it refers to the length scale at which shear is 

4 



comparable to the amplitud of the velocity. Therefore, one expects the formation of clusters 
when A is smaller than the typical interaction diameter, 2R. In other words, the typical 
length scale emerging from the comparison of shear and velocity must be smaller than the 
scale at which we average the velocity of any particle. On the other side, it is clear that when 
Rk. L most of the particles of the system move with the same velocity (all the particles enter 
in the average sum of Q), avoiding the aggregation of the particles. Thus our hypothesis 
for clustering requires that: 

X/2<R<L (7) 

To be specific, in the following the external shear flow is given by v{x) = Uq + 
Vq sin (tux/L), with L the system size (which we take L = 1 so that all length-scales are 
measured in units of L), Uq, Vq positive constants, and uj = 27m, n = 0,1, 2, .... For this flow 



it is not difficult to calculate A = ^/T+'wJJV^ / lo so that on taking ^/l + 2Uq/Vq = 2tt, 
eq. Q becomes 1/n < 2R < 1. 

n 

For a spatial distribution of particles the quantitative characterization of clustering j8|| is 
performed by means of an entropy-like measure 

M 

1=1 

where M is the number of boxes in which we divide the system, and rrii is the number of 
particles in box i. One has that < S]\j < InM, such that Sm = is obtained when all 
the particles are in just one of the boxes, and the InM value is reached when rrii = N/M 
for all i (Poisson distribution of particles), i.e., 5*^ decreases when the clustering increases. 
We define the clustering coefficient as Cm = exp(< Hm >t)/M, where < . >t denotes a 
temporal average at long times, so that when there is no clustering Cm ~ 1. In the left 
panel of fig. (0) we fix = 0.1 (much smaller than the system size L = 1) and plot Cm vs n 
observing that the transition to clustering is obtained for n ^ 5, fitting perfectly ((Tj). In the 
right panel we take n = 10 and plot Cm vs R observing the two transitions indicated in ((7j). 
In figure El we plot the spatial distribution of particles (in the left panel we plot the initial 
distribution) in the regime of clustering at time t = 16 (right) for R = 0.1 and n = 10. Here 
one sees that the particles tend to aggrupate following the sinusoidal fiow. 

Similar results are obtained for other shear fiows. E.g., for the linear shear given by 
v{x) = Fx if X G [0, 1/2] and v{x) = F(l — x) for x G [1/2, 1], A is 1/VT2 so that altering 
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the features of the external flow (the values of F) the aggregation properties of the system 
for fixed R are not changed. However, transitions between non-clustering and clustering 
distributions are observed by varying R. 

In the next section we explain analytically the transition to clustering observed in the 
numerics. This is done by deriving the density evolution equation for the system of particles. 



III. CONTINUUM DESCRIPTION IN TERMS OF THE DENSITY OF PARTI- 
CLES. LINEAR STABILITY ANALYSIS 

A continuum theory can give further insight on the model. The process to obtain 
n 

it is standard f9], and we just present here a sketch: define the particle density as 
p(x, t) = XliliPi(^'^) — 12iLi^i'^ii^) ~ x), then use an arbitrary function /(x) defined 
on the coordinate space, and take the time derivative on both sides of the obvious relation 
/(xj(t)) = / (ix pi(x,t)/(x). Finally, using ijx-x,(t)|<i? '^^ Pi^^t) = ^ni^ one arrives to 



atp(x, t) + V> 



0. (9) 



Note that we have maintained the time dependence of the velocity field to refiect the gen- 
erality of the approach. Eq. Q can be simply read as that the density of particles is driven 
by the effective velocity \'^^^{'x,t) = i|r_x|<_R'^^ v(r,t) p{r,t)/ i|r_x|<i? '^'^ P('^5'^)5 whose de- 
pendence on the density reveals the self-consistent character of the model. Note also the 
two trivial limits: a) i? — > or passive limit, v'^-^-^(x, t) v(x, t), and h) R ^ 1 {L = 1) for 
which (x, t) 1/N J drp{r, t)v(r, t), i.e., the average velocity of the system of particles, 
which is the same for all of them (and constant for a time- independent velocity field). 

Next we make a linear stability analysis of the stationary homogenous solution, pq, of 
eq. 0. We first write p(x, t) = po + ^^(x, t) where e is a small parameter, and ip{'x,t) the 
space-time dependent perturbation, and substitute it in eq. 0- To first order in e, using 
incompressibility of the fiow and denoting . = i|i._x|<_R- obtain 

dtiJ + ^g{^) ■ Vx^ + ^ Vx ■ Jb dv v(r, t)^(r, t) 

-^g(^).J^dT^{r,t) = 0, (10) 

with g(x) = J^dr v(r, t). Though linear, the above expression is still rather complicated 
since it is non-local in space. For the sinusoidal shear fiow (taking for simplicity and without 
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lost of generality Uq = 0), we have that g(x) = y2nR/u! sm{ux)Ji{u!R) with y a unitary 
vector in the y direction, and Ji the first order Bessel function, so that eq. ([nH) becomes 

dt^ + ^^^^^^ sm{ujx)dyij 
_?Y^ sin{ujx)dy[j^ dr ^{r, t)] = 0, (11) 



where r = (r^;, r^). 

We are mainly interested in the clustering transition driven by the relative values of A and 
R, so that we next consider the limit i? << 1. It is very important to note that, to be the 
expansion in R consistent, all length-scales of the system must also be very small compared 
with the system size. Specifically, X = l/u << 1, so that in particular, one cannot expand 
the sm{ujr^) in the integrand in eq. (jllj) . Let us detail the calculations. The two integrals 
appearing in eq. have the approximations (for simplicity of notation we skip the time 
dependence) : 

la = /b dripir^, Ty) = Jj^,|^^ rfrV« + x,r'y + y) 
~ !\v'\<R c^r'[^(x, y) + r'^dyi){x, y) + r'ydyi){x, y)] 

= nR^ij{x,y) + Q{Ry, (12) 
h = /g dr sm{ujr^)tp{r^, Ty) 
= I\r'\<R sin(w<. + ujx)^{r'^ + x, + y) 
~ !w\<R sin(c^< + ujx) [%l){x, y) + r'^dyip + r'^dy^j] 
= HlZii^ sin(cux)^ + ^""^'j^^""^) cos{ujx)d,i; + QiR^). 

(13) 

Here 0(i?^) indicates terms of order R^ and superior. After substituting expressions p2|) - 
(fT^ in eq. the evolution of the perturbation in the small R limit (or better, when the 
typical length scales of the problem are small) is finally given by 



+ sin(.;a;)a,^ + cos(.;x)a^,^ = 0, (14) 



where we have neglected terms of order R. 

Two fundamental features further simplifies the analysis: a) the coefficients are periodic 
in the spatial coordinates so that Floquet theory can be apphed, and b) the coefficients 



are independent of the y coordinates so that plane waves are solutions on the y direction. 
Therefore we make the ansatz: 

oo 

V^(x,y,t)=e^*+*^^+^^^ ^ 0^^e", (15) 

m=— oo 

where, because of periodic boundary conditions, Cj = 2ttp {p = 1,2,...), K = 27rp 
(p = 1,2,...), and 0^ are complex coefficients. K is restricted to the first Brillouin zone 
determined by —uj/2 < K < u/2, and u is not bounded. 

If any of the eigenvalues A is positive then the perturbation grows (the homogenous 
solution is unstable) and clustering emerges in the system. Thus we look for the conditions to 
have A > 0. Using the exponential formula for the sine and cosine functions, and substituting 
expression (fT3j) in eq. ()14|) we obtain, after grouping the terms with the same exponential 
argument, 

Am0m + (pm-i[ai - (3m] + 0m+i[a2 - (3m] = 0, (16) 

with ai = aiCj/2 — a2CjK/2 + a2UjCj/2, = —aiCj/2 — a2CjK/2 — a2UjCj/2, and (i = a2UJU}/2, 
with the notations ai = 2VoJi{ujR) / (uR) and a2 = 4:VoJ2{ujR)/uj. 

For a simple theoretical analysis we just consider the three Fourier modes m = 0, ±1 and 
neglect the rest. Diagonalizing the corresponding 3x3 matrix of coefficients of the system 
in eq. (jlfij) we obtain three eigevanlues, one zero and the other two given by 

A±{K) = ±^JsJlK^R^ - AUJRJ2J1 - 2Jl (17) 
ujR ^ 

where the Bessel functions, Ji and J2, are evaluated at uR. The expression for A+ is 
cuadratic in K with a positive coefficient for the term in K^, so that taking into account 
that —uj/2 < K < uj/2, the inestability is obtained when A+(i^ = uj/2) is positive, i.e., 

J2{ujRfuj'^R^ - 2ujRJ2{ujR)Ji{ujR) - Ji{ujRf > 0. (18) 

Numerically one solves the above inequality and obtains that the condition for instability 
is ujR > 2.5, which, despite the many approximations made to derive it, fits well with the 
numerical result l/u < 2R (eq. ((Tj) for Uq = 0). We have checked that the above result 
is improved by including more modes in the Floquet analysis. In particular, considering 
m = 0, ±1, ±2, the final condition for the maximum exponent to be positive is uR > 
1.32. We believe that in the limit m — * 00 the numerical result is approached. Therefore 
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this analysis confirms that the derived continuum description eq. Q properly describes 
the discrete interacting particles model, and that the clustering emerges as a deterministic 
instability of the density equation. 

IV. SUMMARY 

In this work we have proposed a very simple model for an ensemble of particles self- 
consistently driven by an external shear flow. Despite its simplicity the model shows a very 
interesting behavior where a transition to grouping of particles are observed. An hypothesis 
for the appearence of the clustering has been presented. It esentially says that the clustering 
appears when the length scale that comes from the comparison of the shear flow and the 
velocity field amplitudes is smaller than the typical interaction radius of the particles. This 
hypothesis has been numerically checked and also a continuum description has been derived 
that confirms it. 

A more realistic interaction of the particles, for example decaying with the distance within 
i?, is planned to be studied in the future. Also, it will be interesting a detailed study of the 
role of a noise term in the dynamics of the particles (which in the continuum description is 
a difussion term), and the analysis when a chaotic flow is considered. 
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